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ABSTRACT 

We have measured an extensive set of UBVRIJHK colors for M31 globular clusters (Barmby 
et al. 2000). We compare the predicted simple stellar population colors of three population 
synthesis models to the intrinsic colors of Galactic and M31 globular clusters. The best-fitting 
models fit the cluster colors very well - the weighted mean color offsets arc all < 0.05 mag. The 
most significant offsets between model and data are in the U and B passbands; these are not 
unexpected and are likely due to problems with the spectral libraries used by the models. The 
metal-rich clusters ([Fe/H] > —0.8) are best fit by young (8 Gyr) models, while the metal-poor 
clusters are best fit by older (12-16 Gyr) models. If this range of globular cluster ages is correct, 
it implies that conditions for cluster formation must have existed for a substantial fraction of the 
galaxies' lifetimes. 
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ral zeroth-order test of compatibility between the 
models and reality (Huchra 1996). In this paper 
we compare the broad-band UBVRIJHK colors 
predicted by three modern population synthesis 
models with the colors of Galactic and M31 GCs. 
Our observational database is the first one to in- 
clude extensive coverage of the JHK bandpasses, 
and the first with spectroscopic metallicities for 
all clusters. We use the cluster metallicities to 
bin the clusters for comparison to the appropriate 
models. In this way we determine the cluster-to- 
model offset separately for each color and avoid 
the ambiguity in comparing model and cluster 
colors in two-color diagrams. 

2. Input data and comparison procedure 

For M31 clusters, the observational data are 
from the Barmby et al. (2000) catalog. For 
the Galactic clusters we obtained optical colors, 
metallicity, and reddening from the June 1999 
version of the Harris (1996) catalog, and IR colors 
from Frogel et al. (1980), as reported in Brodie 
& Huchra (1990) (but with the reddening cor- 
rection appHed by Brodie & Huchra (1990) re- 
moved and the reddening values in Harris (1996) 
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Introduction 



Predicting the integrated spectral energy distri- 
butions of stellar populations is important in the 
solution of many problems in astronomy, from de- 
termining the ages of globular clusters to modeling 
counts of faint galaxies at high redshift. Beginning 
with the early work of Tinsley (1968), successive 
generations of modelers have combined the best 
available data on stellar structure and evolution 
to predict the appearance of the combined light of 
generations of stars. Although the subject of pop- 
ulation synthesis has a long history, it is an active 
area of research: synthesis techniques and many 
of the input data (isochrones, opacities, spectral 
libraries) continue to be improved. 

There is good evidence that globular clusters 
(GCs) are internally homogeneous in age and 
metaUicity (Heald et al. 1999; Stetson 1993). GCs 
are the best observational analogs of modelers' 
'simple stellar populations', i.e. populations of 
stars formed over a short time out of gas with 
homogeneous chemical composition. Broadband 
colors are among the simplest predictions of pop- 
ulation synthesis models, so comparing the mod- 
els' predicted colors to cluster colors is the natu- 
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used instead). Wc dcrcddened the clusters' col- 
ors using the values of Eb-v given in the cat- 
alogs and the Cardelli et al. (1989) extinction 
curve for i?y = 3.1. For M31, wc excluded clus- 
ters where the error in the spectroscopic metal- 
licity was ctjfc/h] > 0.5, and clusters suspected 
of being young on the basis of strong Balmer ab- 
sorption or blue B — V colors (see Barmby et al. 
2000). For both galaxies, we excluded clusters 
with Eb^v > 0.5; there are 103 M31 and 85 
Galactic clusters in the final sample. Photometric 
data is not available in all bandpasses for all clus- 
ters: only about two-thirds have measured R and 
I, and less than half have H. 

We compare the cluster colors to those for sim- 
ple stellar populations of ages 8, 12, and 16 Gyr 
from three sets of models: those of Worthey^, 
Bruzual and Chariot (hereafter BC) (both the 
Worthey and BC models are reported in Lei- 
therer et al. 1996), and Kurth et al. (1999) (here- 
after KFF). Although model colors are tabulated 
in smaller age increments (typically 1 Gyr), ini- 
tially it is more reasonable to use the models as a 
rough guide to relative ages rather than attempt- 
ing to derive precise cluster ages from them. The 
Worthey models are computed at [Fe/H] values 
of -2.0,-1.5,-1.0,-0.5, and -0.25 dex, and the 
BC and KFF models are computed at [Fe/H] val- 
ues of -2.33 (KFF models only), -1.63, -0.63, 
and —0.32 dex. We compared clusters to both 
the Salpeter IMF (Worthey's 'vanilla' models) and 
Scalo (1986) (Miller & Scalo 1979, in the Worthey 
models) IMF version of the models. Worthey 
(1994) finds that some of his model colors have 
defects (e.g. B-V is too red by 0.04-0.06™ due 
to problems in the theoretical stellar atmospheres 
and the color-temperature calibration), but the 
sizes of these defects are not well-determined so 
we do not correct for them. Figure 1 shows data 
and models in two frequently-used two-color dia- 
grams. 

Since the models are computed at discrete val- 
ues of [Fe/H], we use the spectroscopic metal- 
licities of the clusters to compare only clusters 
with comparable metallicities (±0.25 dex) to each 
model. The Galactic cluster metallicities given in 

The version we used updates the Worthey (1994) models 
by including a more realistic treatment of the horizontal 
branch for [Fe/H] < -1.0. 



Harris (1996) are on the Zinn & West (1984) (ZW) 
metallicity scale, and the M31 cluster metallicities 
are also tied to this scale through the calibration of 
Brodie & Huchra (1990). Recent work (Carretta 
& Gratton 1997; Rutledge et al. 1997) suggests 
that the ZW scale may be non-linear at both high 
and low metallicities. We retain the ZW scale in 
this paper because we found that using the Car- 
retta & Gratton (1997) scale to assign clusters to 
model comparison bins made little difference in 
our results. We caution, however, that the effect 
of changing the metallicity scale is unknown for 
the [Fe/H] = —0.25 model bin. The transforma- 
tion from the ZW to CG scales is only defined for 
[Fe/H]zw < —0.5, the lower limit of this metallic- 
ity bin. 

We calculated the mean offsets between model 
and cluster colors (referenced to V) for each metal- 
licity bin; Figures 2-3 show some representative 
comparisons. We plot A{X — V) for all bandpasses 
X to make clear the differences in spectral energy 
distributions between models and data; we remind 
the reader that the offsets for bandpasses redward 
of V thus have the opposite sign from the usual 
colors. One general characteristic of the models 
visible in the Figures is that younger- aged models 
predict bluer colors. The exception is the KFF 
Scalo model for [Fe/H] = —1.63, which predicts 
only very small color differences (< 0.0 1'") be- 
tween ages of 12 and 16 Gyr. The effect of the IMF 
on the colors appears to depend on both metallic- 
ity and age. For the Worthey [Fe/H] = —1.50 
models, Miller-Scalo IMF colors are redder than 
Salpeter model colors at all ages, but for the 
[Fe/H] = -0.50 models, the Miller-Scalo IMF col- 
ors are bluer for 8 and 12 Gyr and almost identical 
for 16 Gyr. BC predict almost no color difference 
between the Salpeter and Scalo IMF models of the 
same age and metallicity. 

A striking feature in Figures 2 and 3 is the 
range of discrepancies between models and data. 
For example, the largest difference between the 
Worthey model with parameters (Salpeter IMF, 
[Fe/H] — —1.50, age 16 Gyr) and the mean colors 
of clusters with -1.75 < [Fe/H] < -1.25 is 0.04"» 
inU-V. The same models with [Fe/H] = -0.50 
are well offset from the data at all colors ex- 
cept U-V; the largest offset is 0.23™ inV-K. 
To determine the best-fitting models, we quantify 
the overall goodness-of-fit for each model/cluster 
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metallicity bin pair as: 

The color differences A{X — V)k are weighted by 
where afe are the standard errors in the mean 
colors of objects in the bin. Table 1 gives the A 
and F values for the best fitting models in each 
metallicity bin. 

3. Discussion 

Table 1 shows that the best-fitting models fit 
the data quite well, with typical color offsets of 
0.02 - 0.03™. The two bandpasses with the most 
significant offsets are U and B: the models are 
too blue in U — V and too red in B — V. Neither 
offset shows a clear trend with metallicity. The 
offsets are likely not due to systematics in the pho- 
tometric system or in the extinction curve. While 
problems with the photometric systems might be 
expected in the R and / bands (due to conversion 
between the Johnson and Cousins RI systems), 
both data and models use the well-defined John- 
son UBV system. Problems in the reddening curve 
also seem unlikely for the same reasons. We sus- 
pect that the offsets are more likely due to sys- 
tematic errors in the models. The B — V offset in 
particular is likely due to the flux libraries used. 
Both Worthey (1994) and Lejeune (1997) found 
their model B-V colors to be 0.04 - 0.06™ too 
red compared to empirical solar-metallicity spec- 
tra, even after correcting to the empirical color- 
temperature scale. This suggests a possible prob- 
lem with the stellar atmosphere models of Kurucz 
(1995), upon which both hbraries are based. 

The cause of the offset in C/ — F is not as clear. 
This offset is actually worse than it appears: since 
we compute the U — V colors for the BC and KFF 
models as (U-B) + (B-V), the red B-V colors 
compensate for some of the U — B defect, which is 
actually larger than the defect in U — V. Worthey 
(1994) - whose models give U — V directly - finds 
that his model U — V is too blue compared to 
solar neighborhood stars and elliptical galaxies. 
Worthey cites problems with the U fluxes from 
the stellar libraries as a possible cause: model- 
ing the many blended atomic and molecular lines 
blueward of B is difficult, and many of the nec- 
essary opacities are not well determined. This 



cannot be the only cause of model problems in 
U, since the BC and KFF models, which use the 
same stellar library, predict different U — V col- 
ors. The treatment of the horizontal branch in 
the models is another possible source of problems 
in the U — V colors because the HB emits most of 
the blue light. However, systematic problems with 
the model HB color (which depends on metallic- 
ity), would presumably produce a. U — V offset de- 
pendent on metallicity, which we do not observe. 
Observational error is another possible contributor 
to the U — V offset, as many of the U — B colors of 
the M31 clusters are poorly determined (sec Table 
3 of Barmby et al. 2000). Understanding the rest- 
frame U flux of stellar populations becomes in- 
creasingly important when studying high-redshift 
galaxies and global star formation history, and fur- 
ther investigation of the models in this bandpass 
is clearly warranted. 

A secondary result in Table 1 is that age de- 
termines which model best flts the data. Higher- 
metallicity cluster colors are best fit by 8 Gyr 
models, regardless of IMF. Lower-metallicity clus- 
ter colors ([Fe/H]bi„ < -1.00) are best flt by 12 
or 16 Gyr models. The best-fit age depends on 
the IMF for several of the models, but not in 
any systematic fashion. This result is consistent 
with the determinations of relative ages for Galac- 
tic clusters by Rosenberg et al. (1999). These 
authors determined relative ages of 35 Galac- 
tic globular clusters from a homogeneous set of 
V, V — I color-magnitude diagrams. They com- 
pared theoretical isochrones with the observa- 
tional CMDs to determine ages using two indepen- 
dent methods. They found that the clusters with 
[Fe/H]cG > —0.9 were ~ 17% younger than clus- 
ters with [Fe/H]cG < ^1-2, with the intermediate- 
metallicity clusters showing a ~ 25% age disper- 
sion. These results are model-dependent, as are 
ours, but the results' similarity implies that ei- 
ther there is a real difference between metal-rich 
and metal-poor clusters or there is a systematic 
problem in the models in one of the metallicity 
regimes. 

What possible systematic errors in our input 
data or comparison procedure could produce the 
result that the metal-rich clusters are younger? 
We redid the comparison procedure considering 
the clusters of each galaxy separately, and still 
found younger ages for the most metal-rich clus- 
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ters. Although M31 has a greater proportion of 
the metal-rich clusters, younger ages are found for 
both M31 and Galactic metal-rich clusters. Co- 
hen & Matthews (1994) suggest that the spec- 
troscopic metallicities of the most metal-rich M31 
clusters measured by Huchra et al. (1991) are sys- 
tematically too high. If this is true, the clusters 
would appear too blue compared to old, higher- 
metallicity models and the best-fit model would 
be younger. We compared the metal-rich M31 
clusters to the Worthey [Fe/H] = —1.0 models, 
and the best-fitting model had age 16 Gyr. How- 
ever, the goodness-of-fit was better for the young, 
metal-rich models than for the older, more metal- 
poor model, so we conclude that younger ages 
arc still favored for these clusters. Overestimat- 
ing the reddening of the metal-rich clusters would 
make the derived intrinsic colors too blue and yield 
younger ages. This seems unlikely, given that the 
color-metallicity relations for Galactic and M31 
clusters match well throughout their metallicity 
range (sec Barmby et al. 2000), and the methods 
of reddening determination for M31 and Galactic 
clusters are different. 

If the detection of younger ages for metal-rich 
globular clusters is real, it has implications for 
galaxy formation. A range of GC ages implies 
that GC formation took place over an extended 
period of time. Conditions for GC formation were 
not particular to the early universe, an assertion 
supported by observations of 'proto-globular' clus- 
ters in present-day merging galaxies (e.g. Zepf 
et al. 1999). More precise knowledge of the distri- 
bution of cluster ages in each galaxy would be ex- 
tremely useful in understanding cluster system for- 
mation. If the age distribution is continuous, the 
relation between age and metallicity might hold 
clues as to what factors controlled the cluster for- 
mation rate. If the age distribution is bimodal - 
with most clusters old and coeval and the remain- 
der younger and coeval - then some event must 
have triggered the second episode of GC forma- 
tion. Perhaps the younger clusters were stripped 
from or accreted along with satellite galaxies of 
M31 and the Galaxy. 

4. Conclusions 

Comparison of three sets of population syn- 
thesis models with integrated colors of M31 and 



Galactic globular clusters shows that the mod- 
els reproduce the redder average cluster colors to 
within the observational uncertainties. The poorer 
agreement in U — V and B — V is likely due to sys- 
tematic errors in the models. Younger-age mod- 
els are required to best match the colors of the 
metal-rich clusters, consistent with the findings of 
Rosenberg et al. (1999) that the most metal-rich 
Galactic clusters are younger than the bulk of the 
globular cluster population. A range of ages for 
globular clusters implies that conditions for clus- 
ter formation were not restricted to the early uni- 
verse. The cluster age distribution has important 
implications for galaxy and globular cluster sys- 
tem formation, and attempts to determine it more 
precisely are needed. 
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Fig. 1.— {B - V)o vs. {U - V)o and {V - K)o 
for M31 globular clusters (triangles) and Galac- 
tic GCs (sqtiarcs). Lines are population synthesis 
models of ages 8 Gyr (bluer colors) and 16 Gyr 
(redder colors): Worthey (solid), BC (dashed), 
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Fig. 2. — Color offsets A{X — V) (data— models) for Salpeter IMF. {AV = is plotted to emphasize that the 
models arc normalized to the data at V.) Solid lines: 16 Gyr models, dotted lines: 12 Gyr models, dashed 
lines: 8 Gyr models. The —1.5 and —0.5 bins in [Fe/H] are Worthey models; the —1.63 and —0.63 [Fe/H] 
bins are BC models. Error bars (plotted only on the 8 Gyr models for clarity) are the standard errors of the 
mean cluster colors and do not include observational uncertainties. 
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Fig. 3.— Color offsets A{X - V) or A{V - X) (data-models) for Scalo or Miller-Scalo IMF. Symbols same 
as Figure 2. The —1.5 and —0.5 bins in [Fe/H] arc Worthey models; the —1.63 and —0.63 [Fe/H] bins are 
KFF models (which have no A(J — V), A{H — V) since they do not predict these colors). 



Table 1 

Color offsets for best-fitting models 



[Fe/H] 


model 


F 


Ngc 


A(!7 - V) 


A(B - V) 


A(y - R) 


A{V-I) 


A(V-J) 


A{V - H) 


A(y - K) 


-0.25 


W8sc 


0.023 


15 


23.4 ±2.4 


-0.4 ±0.9 


-3.0 ±1.1 


-0.3 ± 1.0 


3.2 ±2.6 


-1.4 ±6.0 


-6.0 ±4.9 


-0.32 


K8sc 


0.022 


17 


11.5 ±2.2 


-2.5 ±0.8 


-2.1 ± 1.0 


0.2 ±0.9 






-3.9 ±4.6 


-0.50 


W8sc 


0.013 


30 


22.8 ±3.1 


-0.3 ±0.8 


-1.0 ±0.7 


1.2 ±1.3 


4.6 ±4.6 


2.8 ±5.2 


-3.3 ± 7.2 


-0.63 


K8sp 


0.012 


35 


12.2 ±3.2 


-1.0 ±0.9 


-0.6 ±0.6 


2.4 ± 1.2 






2.1 ±5.9 


-1.00 


W12sc 


0.014 


48 


13.7 ±2.2 


-2.3 ±1.0 


0.0 ±0.5 


-1.1 ± 1.2 


3.2 ± 1.8 


6.9 ±2.5 


4.3 ±2.1 


-1.50 


W16sp 


0.016 


75 


4.1 ± 1.4 


-3.0 ±0.6 


-0.3 ±0.5 


-1.0 ±0.8 


2.0 ± 1.7 


-0.7 ±3.5 


2.0 ±2.3 


-1.63 


K16sp 


0.019 


77 


5.1 ± 1.2 


-0.7 ±0.5 


1.7 ±0.4 


2.8 ±0.7 






10.8 ±2.1 


-2.00 


W16sp 


0.027 


37 


8.2 ± 1.2 


-2.6 ±0.6 


-0.3 ±0.8 


-1.6± 1.2 


-6.7 ±2.6 


-5.2 ±3.7 


-5.8 ±3.5 


-2.23 


K16sp 


0.027 


12 


4.1 ± 2.3 


1.8 ± 1.3 


2.3 ±2.1 


3.4 ±2.0 






-7.1 ±6.2 



Note. — Units are hundredths of a magnitude. The capital letter in the model column indicates the best-fitting model (Worthey, 
EC or KFF), the number is the model age in Gyr, and 'sp' or 'sc' indicate Salpeter or Scalo IMF, respectively. 
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